N93-80847 


The Initiation, Propagation, and Effect of Matrix Microcracks 
in Cross- Ply and Related Laminates * 

John A. Nairn, Shoufeng Hu, Siulie Liu, Jong Bark 

Material Science and Engineering Department / / 

University of Utah, Salt Lake City , Utah 84112 



ABSTRACT 

Recently, a variational mechanics approach has been used to determine the thermoelastic 
stress state in cracked, [0„/90 m ], laminates (J. A. Nairn, J . Comp . Mat., 23, 1106 (1989)). This 
paper describes a generalization of the variational mechanics techniques to handle other cross- 
ply laminates ([90 n /0 m ] a ), related laminates ([±#/90 m ] a , [90 m / ± B\ s , etc.), and to account for 
delaminations emanating from microcrack tips. Microcracking experiments on Hercules 3501-6/AS4 
carbon fiber/epoxy [90 n /0 m ] a laminates show a staggered cracking pattern. These results can be 
explained by the variational mechanics analysis. The analysis of delaminations emanating from 
microcrack tips has resulted in predictions about the structural and material variables controlling 
competition between microcracking and delamination failure modes. 

INTRODUCTION 

Many observations have confirmed that the initiation of damage in multidirectional laminates is often 
by microcracks in the off-axis plies that run parallel to the fibers in those plies [1 — 8] . These microcracks have 
typically been studied in cross-ply laminates in which the cracks form in the 90° plies [1-8]. Microcracks 
form during static testing [1-8], during fatigue testing [3,9,10], and during thermal cycling [11]. Because 
microcracks cause a reduction in stiffness [3], a change in the thermal expansion coefficient [12,13], and 
provide sites for the initiation of delaminations, it is important to gain a quantitative understanding of the 
formation and propagation of microcracks. Two important factors that must be understood are 

1. The formation of microcracks in 90° plies is dependent on laminate structure. That is, the cracking 
process depends on the thickness of the 90° layers, the support provided to the 90° layers by other plies, 
and whether the 90° plies are on the inside or are adjacent to a free edge. 

2. The residual thermal stresses in the 90° plies of cross-ply laminates are typically tensile and of sufficient 
magnitude to influence the microcracking process. 

STRESS ANALYSIS TECHNIQUES 

Several attempts have been made at the stress analysis of cross-ply laminates using simplistic analyses 
(e.g. shear-lag analyses [2,4,6,7,14,15]). While these simplistic models often yield qualitatively reasonable 
results, they are unsuitable for a thorough understanding of microcracking. Hashin [16,17] was the first 
to apply variational mechanics techniques. He solved for the modulus [16,17] and the thermal expansion 
coefficient [18] of [0 n /90 m ] a laminates as a function of microcrack density. Nairn et al. [19,20] have extended 
Hashin s analysis to include residual thermal stresses and to the fracture mechanics analysis of microcracking 
in [0„/90 m]j laminates. This paper describes the variational mechanics approach, how it can be extended to 
more difficult problems involving [90„/0 m ] a , [±0/9O n ] a , and [9O n /±0] a laminates, and how it can analyze 
delaminations emanating from microcrack tips. In this section, we outline the variational mechanics stress 
analysis techniques. 

Consider a laminate plate with the x axis parallel to the zero degree fibers, the y axis parallel to the 
90 fibers, and the z axis normal to the plane of the plate. When the 90° plies have through-the-width 
microcracks (as is invariably observed in static testing), a two-dimensional analysis of the x-z plane suffices. 
The x-z plane, or laminate edge, can conveniently be divided into multiple layers (n layers). The layers may 
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be assigned to individual plies, ply groups, or portions of a ply, depending on the nature of the problem 
being solved and on the desired accuracy. For the stress state within any layer, we make one and only one 
assumption — that the x-axis tensile stresses depend only on x and are independent of z. A general stress 
state for the i th layer that fulfills this assumption and equilibrium is [19,20] 

4*2 = <4o - MO W 

49 = Ai/<«, 0 < 2 ) 

<*? = -*?*(*. 0 < 3 ) 


where superscript (i) denotes stress in the i th layer, <r§ is the x axis tensile in the uncracked laminate, 
= ti/t,, ipi(£) is a function of the dimensionless x direction coordinate = x/<»J, fi and are functions 

of £ and of the dimensionless z direction coordinate (c = ( z ~ 4 0 )/**')’ and 4 0 ls the z coordlnate at the 
start of the i tk layer. In these equations, f, is the thickness of the i ,h layer and t, is the thickness of any 
arbitrarily chosen normalization thickness. The functions /,, gi, and V« are interrelated by 

dfj _ dgj 

di ~ dC 
„ _ H'i _ dfi 
^ ~ di d( 

We consider a unit cell of damage of a multilayered sample that extends from x = a to x = -a or from 
£ _ p — aft, to ^ — —p. We further consider each layer to be orthotropic with at least one symmetry axis 
aligned with either the x or 2 axes. Generalizing the two-layer results from Ref. [19] to n layers, the total 
complementary energy in the unit cell of damage per unit depth (y direction dimension) can be written as 


(4) 

(5) 



where E xi E z i, Gi, a xi , and a zi are the mechanical properties of the i th layer, being, respectively, x and 
z direction tensile moduli {E„ and E„), in-plane shear modulus (G„), in-plane Poisson’s ratio K), and 
x and z direction thermal expansion coefficients. The term AT = T 0 - T, is the temperature differential 
between the stress-free temperature (T 0 ) and the specimen temperature (T,). The term r 0 is a constant 
energy term that does not enter energy minimization procedures. 

The analysis procedure for a large variety of problems are similar: 

1. From observation of experimental results, model microcracking damage by a unit cell of damage extend- 
ing from x =■ — a to x = a. 

2. Divide the laminate into n layers where the divisions are chosen for sufficient accuracy. 

3. Using boundary conditions and stress continuity conditions, express the fi and ^functions in terms of 
the ifti functions and explicit functions of £ (this step is always possible). 

4. Rewrite the total complementary energy in integral form involving only the V>» functions. 

5. Using the calculus of variations, minimize the complementary energy to solve for the ipi functions. The 
principals of variational mechanics state that these functions will provide the best approximation to the 

true stress state. 

6 With knowledge of the Vi functions, it is possible to find the sample modulus [16,17,19], the thermal 
expansion coefficient [18,19], the total strain energy [19,20], and the energy release rate due to the 
propagation of damage [19,20]. 
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When n is large, solving for all the V» functions will necessarily involve numerical calculations. The solution 
of such problems can be viewed as a complementary energy based finite element (or layer) analysis. When 
n is small (n < 4), however, it is often possible to eliminate all but one or two of the Vi functions. In 
such situations, the resulting calculus of variations problem can be solved in closed form. Fortunately, 
many interesting and relevant cross-ply laminate problems can be expressed with four or fewer layers. The 
remainder of this paper outlines the solutions to some of these problems. 


MICROCRACKING IN [0„/90 m ], LAMINATES 

When loaded in tension, [0 n /90 m ], laminates fail in a nearly periodic array of through-the-width mi- 
crocracks in the 90° plies. The unit cell of damage is shown in Fig. 1A. Because the problem is symmetric 
about the midplane, the laminate can be divided into two layers of thicknesses t x and t 2 (see Fig. 1A). Layer 
1 is half the 90° ply group and layer 2 is one of the 0° ply groups. We need to solve for Vi and V2, one of 
which can be eliminated by force balance. As shown in Ref. [19], the total complementary energy in the unit 
cell of damage per unit depth can be written in terms of Vi as 


r — r 0 + 1 


1 f [C-iV? + C2 V 1 V \ + C3V” 2 + C 4 V[ 2 — 2Aa A'/' t/’i ] dS, 

J — p 


where Aa = 07- — a A , and 


(7) 
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Xva 
3 E a 


( 8 ) 


Here, E, G, v and or denote tensile and shear moduli, Poisson’s ratio, and thermal expansion coefficient, 
subscripts A and T denote axial and transverse properties of the ply material, and X = f 2 /<i. 

Minimization of the complementary energy in Eq. (7) has been accomplished in Ref. [19], We quote 
some useful results. Consider a sample with N unit cells of damage characterized by crack spacings pi, p 2 ,. 
Pat- The sample compliance is [16,19] 


C — Cq + 


KiCsLE* eL xM 
B 2 WE e 2 EL Pi 


(9) 


where W is the sample width (y-direction dimension), L is the sample length (re-direction dimension), E c is 
the modulus of the uncracked sample, and Co = L/BE C W is the compliance of the uncracked sample. The 
new function y(p) has a physical interpretation as being proportional to the excess strain energy caused by 
the presence of the microcracks. Defining p = (C 2 - C 4 )/C 3 , q = C 1 /C 3 , 


where for Aq/p 2 > 1 


and for 4 q/p 2 < 1 


X(p) = 2 ap (a 2 + /? 2 ) cosh 2op- cos 2/?p 
' 0 cosh ‘lap + a sin 20 p 

xW = - » i )- t, " h / >>tanh ' > f 


0 tanh 0 p — a tanh ap 


for t| > , 

P 2 

for < 1 
P 2 



( 10 ) 

(11) 


( 12 ) 


(13) 


We note that the expression for sample compliance (Eq. (9)) has no adjustable parameters and is applicable 
to any distribution of crack spacings (not just periodic arrays of cracks). Hashin [16] has shown that this 
expression is in excellent agreement with experimental data. 
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Table I: The critical microcracking fracture toughness, G mc , for five different carbon fiber composite material 
systems. Details of experiments given in Ref. [20]. 


Prepreg Material 

Gmc (J/m 2 ) 

Hercules 3501-6/AS4 
Fiberite 934/T300 
DuPont Avimid® K Polymer/IM6 
Fiberite 977-2/T300 
ICI PEEK/AS4 

240 

690 

960 

1800-2400 

3000 


The total strain energy in the N crack intervals can be written as [19] 


U = 


/(Tg t 1 Aa 2 AT 2 \ 
V2 E c BC\ ) 


BWL + (C - Co) 


B 2 W 2 E 2 C 

2E$ 


El , Aa 2 AT 2 \ 

Ef^ C\ ) 


(14) 


where cro is the total stress applied 
cracked sample is [19,20] 


to the laminate. The longitudinal thermal expansion coefficient of the 


(*l = <*L ~ 


C — C 0 Aa 
Co CiE t 


(15) 


where a° L is the longitudinal thermal expansion coefficient of the uncracked sample. Finally, combining 
all expressions, we can derive an expression for the energy release rate due to the formation of additional 
microcracks [19,20] 

where (x(p)) is the average value of x{p) over the N crack spacings and D = ^ is the crack density. We 
note that Eq. (16) is slightly different than the result from Ref. [19] and corrects an error in that paper. The 
correct result is derived in Ref. [20]. 

If we assume that the formation of microcracks is governed by a critical energy release rate or a mi- 
crocracking fracture toughness, G mc , we can use Eq. (16) to predict the initiation and the increase in crack 
density in [0„/90 m ] s laminates. To use Eq. (16), we must evaluate the density derivative of D{ X (p)}- The 
formation of a new microcrack between two existing microcracks is shown iu Fig. IB; the crack is shown to 
form in the middle because that is the location of the maximum stress and because of the tendency of these 
laminates towards periodic crack spacings. By a discrete evaluation of the required derivative 


G m = 2 C 3 ti(2x(p/2)-x(p)) ( 17 ) 

Given a value of the fracture toughness of a material system, Gmct we can solve Eq. (17) for applied 
stress, <To, and predict the crack density as a function of applied stress. This approach was applied in Ref. [20] 
to five material systems; for each material system, two to five cross-ply laminates were used. Some typical 
results for Hercules 3501-6/AS4 carbon fiber/epoxy laminates are given in Fig. 2. For this material, G me 
was found to be 240 J/m 2 . An important result for this and other materials systems is that a single value 
of G mc for a given material system suffices to predict the results from all cross-ply laminates tested. Some 
measured values of G mc are given in Table I. 

The variational mechanics analysis of microcracked [0„/90 m ], laminates has been useful. Without any 
adjustable parameters it gives an excellent prediction of the sample modulus as a function of crack density 
[16]. When implemented into a fracture mechanics analysis of microcracking, it can predict the crack density 
as a function of applied load [19,20], Combining experimental results and the energy release rate expression 
in Eq. (17), it is possible to measure G mc ■ G mc appears to be a useful material parameter characterizing 
microcracking or intralaminar fracture toughness of composite materials. 
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MICROCRACKING IN [90 n /0 m ] a LAMINATES 

The damage process in [90 n /0 m ] a laminates is more complicated than in [0 n /90 m ] a laminates. When 
[90 n /0 m ] a laminates are loaded in tension, the 90° plies on either side develop nearly periodic arrays of 
through-the-width microcracks. Comparing the crack patterns on either side, however, we observe that the 
cracks on one side are shifted by half a crack spacing from the cracks on the opposite side. Thus, any 
given microcrack is located approximately midway between two microcracks on the opposing surface. The 
unit cell of damage for such “staggered” microcracks is shown in Fig. 3A. Because the problem is no longer 
symmetric about the midplane, we must analyze the entire laminate; we thus divide the laminate into four 
layers according the ply groups (see Fig. 3A). Layers 1 and 4 are 90° ply groups and layers 2 and 3 divide 
the 0° ply group into two layers. We need to solve for to ^ 4 , two of which can be eliminated by force 
balance and transverse stress free-edge effects. Finding the complementary energy for this system is a much 
more difficult than for the [0 n /90 m ] a laminates. The details will be described in a future publication. In this 
paper we give the key results. The total complementary energy in the unit cell of damage per unit depth in 
terms of and rp 4 is 


r = r 0 + t\J (c^ 2 + c 2a xx" + c 3a x " 2 + c 4 x 12 - 4A<*atx 

+ C\Y 2 + C^YY" + C 3 Y " 2 + ClY' 2 )dt 

where X = rf>i + xp 4 , Y = ipi — Vm, Cj and C 4 are given in Eq. (8), and the new constants are 


r - ** 4 . 

1 (1+2A) 2 

1 Et + A 3 E a 




C?,a = 


+ ^ (8A 2 + 20A + 15) 


10 Et OOEt 


(7* _ _ _f*L_ i 
2 3 E t + 


va 

E a . 


(1 + 2A)(2 + A) 


3A 


c -2ob + 5ok(^ + 7i + 8 ) c < - = 357 + i W ! 


(18) 


(19) 


Minimizing the complementary energy in Eq. (18) would be an intractable problem if it were not for 
the following symmetry relation between V>i and 



Vh(£ -p) 

Mt + p) 


for £ > 0 
for £ < 0 


(20) 


With this relation, it is possible to minimize the complementary energy in closed form. The details (which 
are more complicated than the [0 n /90 m ] a problem) will be described in a future publication. 

Because the stress analysis of [90 n /0 m ] a laminates is a new solution, we begin with a brief discussion 
of the resulting stresses. The tensile stress <ri l J } or the tensile stress in the 90° ply on the left of Fig. 3A is 
plotted in Fig. 4; the plot is for a Hercules 3501-6/AS4 carbon fiber/epoxy [902/0]* laminate with a crack 
spacing characterized by p = 3. At the two crack faces cr^ = 0 as required by boundary conditions. Midway 
between the two microcracks and directly opposite the crack in the 90° ply on the opposing surface (see 
Fig. 3A) there is a local minimum in tensile stress. This local minimum is caused by a bending effect that 
results from the asymmetric arrangement of microcracks. Two local maxima in stresses are located near 
p = ±1 or at positions roughly ^ and | of the way from the bottom microcrack to the top microcrack. The 
form of the stresses in Fig. 4 can be used to explain the tendency towards staggered microcracks. In Fig. 3B, 
we show new microcracks formed at all local stress maxima. It is observed that the new microcrack pattern 
is equivalent to three unit cells of damage and thus the distribution of the stresses naturally leads to the 
propagation of staggered microcracks. 

We next quote some useful results of the new [90 n /0 m ] a laminate analysis. It is possible the cast these 
results in a form that is reminiscent of the [0 n /90 m ] a laminate results. Consider a sample with N unit cells 
of damage characterized by crack spacings pi,p2>- • • , Pn> The sample compliance is 


2 hCaoLEl EL XaiPi) 

B 2 WEl £* 


t-T 
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where the new function Xa(/>) has a physical interpretation as being proportional to the excess strain energy 
caused by the presence of the microcracks; it is defined by 


Xa(p) 


ld- 


Mil 


(22) 


In Eq. (22), x(p) is defined by Eq. (11), except that C 2a and C 3a replace C 2 and C 3 , and x’(p) is a new 
function. Defining p* = (C 2 — C^/Q, q * = C\ jC 3 , 


X» = 2a*/T (a* 
X'(P) = «*/?* (/?* 2 


2 a* 2 \ cosh 2 a*p 4 - cos 2/?*p 

/ /?* sinh 2 a* p — a * sin 2/3* p 

_ a »2\ 1 

/ /?* tanh a* p — a* tanh (3* p 


4o* 

for -nr > 1 
P* 2 


for 


v 

p* 2 


< 1 


where for 4 q* /p * 2 > 1 

«* = and (3* = ^2^ + p* 

and for 4q* /p* 2 < 1 



(23) 

(24) 


(25) 


(26) 


As with [0 n /90 m ], laminates, the expression for the compliance of [90^/0™], laminates (Eq. (21)) has no 
adjustable parameters and is applicable to any distribution of crack spacings (not just periodic arrays of 
cracks). We do not yet have experimental data that can be used to compare predictions to observations. 

The total strain energy in the N crack intervals and the longitudinal thermal expansion coefficients are 
given by expressions identical to the [0 n /90 m ]^ laminates (see Eqs. (14) and (15)). A difference between 
[90 n /0 m ] 3 and [0 n /90 m ] s laminates occurs because the expressions for compliance, C, that must be used 
with Eqs. (14) and (15), differ. Finally, combining all expressions, we can derive an expression for the energy 
release rate due to the formation of additional microcracks 


Gm = (fr*o- ^^) 2 C 3a < 1 ^(D(x„(p))) (27) 

where (xa(p)) is the average value of Xa(p) over the N unit cells of damage. 

If we assume that the formation of microcracks is governed by a critical energy release rate or micro- 
cracking fracture toughness, G mc , we can use Eq. (27) to predict the initiation and the increase in crack 
density in [90 n /0 m ] 5 laminates. To use Eq. (27), we must evaluate the density derivative of D(x a (p))- If we 
consider the cracking process in Fig. 3B, which perpetuates the staggered arrangement of microcracks, we 
can evaluate the derivative by the methods used for [0 n /90 m ], laminates [19,20]. The result is 


G m = 1 -(^a 0 - ^P) 2 C 3a <i(3 Xa (p/3) - Xa (p)) 


(28) 


Given a value of the fracture toughness of a material system, G mc , we can solve Eq. (28) for applied 
stress, <r 0 , and predict the crack density as a function of applied stress. Our first experiments have been 
with Hercules 3501-6/AS4 carbon fiber/epoxy laminates; the results for a [902/0/902 ]t laminate are given 
in Fig. 5. To be useful, G mc should be a material constant that is independent of stacking sequence. We 
thus fit the [902/0/902 ]t laminate results using the G mc of 240 J/m 2 found during experiments on [0 n /90 m ], 
laminates [20]. The results in Fig. 5 show that Eq. (28) fits the experimental results at high crack density 
and at zero crack density (microcrack initiation) but shows some deviations at crack densities between 0.0 
and 0.4 cracks/mm. The cause of the deviations will require more analysis. We speculate, however, that it is 
related to the discrete evaluation of the density derivative of D(x a (p)} inherent in Eq. (28). When the crack 
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density is low and the microcracks are far apart, the local stress maxima located at ±pf 3 become diffuse and 
move away from ±p/ 3 towards ±p/ 2. In other words, it is only at high crack densities that the distribution of 
the stresses causes cracks to form at ±p/3, that damage is driven towards a staggered cracking pattern, and 
that Eq. (28) gives the true energy release rate. At low crack densities, the cracks will be more random and 
the expression for energy release rate needs to be modified. Experiments which measure the crack patterns 
at low crack densities can be used to confirm this speculation and to suggest ways to modify Eq. (28). 

Although much more difficult than previous analyses, a closed-form variational mechanics analysis of 
[90„/0 m ], laminates has been obtained. Literature results concentrate of [On/flOmJj laminates to the near 
exclusion of [90„/0 m ], laminates. Clearly more experimental results for [90„/0 m ] 3 are desirable. It is 
our expectation that the variational mechanics analysis will provide the tools necessary for a quantitative 
understanding of those results. 

MICROCRACKING IN [±0/9O m ], AND IN [9O„/±0], LAMINATES 

The results in the previous two sections can readily be adapted to handle problems of microcracking in 
90° plies supported by any orthotropic sublaminate. One simple example is microcracking in [±0/9O m ], and 
in [90„/ ±9] t . In these laminates, 6 should be relatively small ( e.g . 0 < 45) or else microcracking of the 90° 
plies may not be the dominant failure mode. The adaptation of the variational mechanics analysis consists 
merely of adjusting the constants C x to C 4 , C 2a , C 3a , and C{ to C\. We note that each of these constants, 
defined in Eqs. (8) and (19), is the sum of two terms. In ail cases, the first term results from the 90° plies 
and the second term results from the 0° plies. To construct a solution for laminates in which the 0° plies are 
replaced by a [±0] sublaminate, we merely replace the mechanical properties of the 0° plies in each second 
term with those of the [±0] sublaminate. 

DELAMINATIONS EMANATING FROM THE TIPS OF MICROCRACKS 

All variational mechanics analyses of cross-ply and related laminates can be extended to account for 
delaminations emanating from the tips of existing microcracks. We illustrate the delamination analysis using 
[0 n /90 m ], laminates. Figure 6 shows the region between two microcracks separated by a distance 2 a having 
delaminations of length d\ and d >2 emanating from the top and bottom microcracks and propagating into 
the region between the two microcracks. 

We consider the right half of the symmetric laminate in Fig. 6. We split the area between the two 
existing microcracks into three regions. Region I is the region within the top delamination or the region 
from x — a — to x — a. Region II is the region between the tips of the delaminations or the region from 
x = ~ a T ^2 to x = a — d\ . Region III is the region within the bottom delamination or the region from 
x = -a to x = -a + d 2 . Within each region we divide the right half of the laminate into two layers. Layer 
1 is the 90° ply group and layer 2 in the 0° ply group. 

We begin with the stresses in the 90° plies in regions I and III. By symmetry, the shear stress at 
the sample middle must be zero (cr^(0) = 0). Because of the requirement for stress free surfaces on the 
delamination crack face, the shear stress and the normal stress at the right edge of the 90° plies must also 
be zero (cr^(l) = 0 and cr^(l) = 0). From these boundary conditions and the form of the stresses in 

Eqs. (l)-(3) which assume that the normal stress is a function of x only, the stresses in regions I and 
III are uniquely determined to be 


= 0 

4 ^ = <7 (1) - 


= 0 


1 + A 


4 * = 


<t(V = <r( 2 ) = 

U XZ U 2Z 


(29) 


The boundary conditions in region II are unaffected by the length or even the presence of the delaminations. 
With or without the delaminations, the x direction tensile stress and the shear stress are both zero at the 
top and bottom of region II. The stresses in region II are therefore identical to the stresses calculated in 
Ref. [19] except that the crack spacing needs to be adjusted. We introduce the dimensionless quantity 


di -f d 2 

6= — t r 


(30) 
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The stresses in region II are identical to the stresses calculated in Ref. [19] for two microcracks whose spacing 
is characterized by p — 6. 

With the stresses completely determined, we can follow procedures similar to those in Ref. [19] and 
evaluate properties of microcracked and delaminated cross-ply laminates. Consider a sample with N unit 
cells of damage characterized by crack spacings pi,P2>- • ■ > Pn an d by extents of delamination within each 
crack interval of 61,62 >■ ■ ■ Av- The sample compliance is 


_ L — d f 2t 1 C 3 LE*'Ex(pi-6i) \ , ± + * L \ 

" L V B'WEl £(pi - 6i) J L \ X BE a W) 


(31) 


where d = * s total sample delamination length. Although we have no experimental evidence 

to verify Eq. (31), its limiting values are appropriate. In the limit of no delaminations, or d — 0, Eq. (31) 
reduces to Eq. (9). In the limit of complete delamination, or d = L, the sample compliance reduces to 


Coo — 


1 + x 


X be a w 


(32) 


This is the compliance for a sample in which the 90° plies carry no load. In terms of compliance, the total 
strain energy is found to be 


" - (m- * (L ~ d)BW - + ikrr dBW 

+ (c — Co — -j- (Coo — Co)) 
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and the thermal expansion coefficient is found to be 

n d Ac* Ac* 


<*L = <*L ~ 


L A £4 Ci C x C q Et 


(c-Co-| (Coo - Co)) 


(34) 


The above two-dimensional analysis can immediately be used to find the energy release rate for the 
growth of a through-the-thickness delamination — CY The energy release rate is 


r - du 

Gd ~~dA 


const. disp. 


1 du 

'2 W dd 


(35) 


const, disp . 


where A is total delamination area which is equal to 2 dW . We differentiate Eq. (33) realizing that <To and 
C will depend on d. The result after much simplification (details can be found in Ref. [20]) is 


Ga 


Ej- 

~E, 


-<T 0 


AaATY 
Ci ) 


C 2 t x Y D (pJ) 


(36) 


where the function Yo(p,S) depends on p = (pi,pi, ■ ■ ■ ,Pn) and on 6 

m^) 


= (6i ,6 2 , • • • , ^jv)- R is given explicitly 


d Ex(p. ~ h) 
dd E(Pi “ 6 >) 


(37) 


Imagine a sample with N crack intervals characterized by dimensionless crack spacings and delamination 
lengths defined by the vectors p and 6. Imagine further some delamination growth confined to the k th crack 
interval. The Y D (p, 6) function for this delamination growth simplifies to [20] 


W.«) = ^(x , (0)-x , (pfc-«*)) 

The differentiation of \(p) is with respect to the variable p. 


(38) 
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If we assume that delamination will propagate when G d > Gdc where Gdc is a critical energy release 
rate or interlaminar fracture toughness, we can use Eq. (36) to predict the initiation and growth of micro- 
crack induced through-the-width delaminations. One possible problem of this approach is that it uses total 
energy release rate instead of partitioning it into the various modes — mode I, mode II, and mode III. We 
suggests that the assumed ^-direction propagation of a through-the-width delamination will be mostly mode 
II deformation. The distinction between using total energy release rate and the mode II part of the energy 
release rate will therefore be minor. Experimental results are required to determine the adequacy of using 
total energy release rate as a failure criterion. For the purpose of this discussion, we assume that the total 
energy release rate failure criterion is valid. 

We can use the delamination analysis to discuss competition between microcracking and delamination. 
The first form of damage in cross-ply laminates is always microcracking. Once the first microcrack has 
formed, we can ask if a delamination will initiate at that microcrack or if instead another microcrack 
will form. Comparing the energy release rate for microcracking (Eq. (17)) to the energy release rate for 
delamination (Eq. (36)), the predicted failure mode will depend on the values of G mc and G dc and on the 
relative magnitudes of 2\(p/2) - \{p) and ^ (x'(0) - *'(/?)). Assuming G mc and G dc are the same (they 
both represent crack growth through the matrix, albeit possibly by different fracture modes — mode I vs. 
mode II), the predicted failure mode can be determined by plotting the latter two quantities. Figure 7 plots 

2x(p/2) “ x{p) an d ^ (x 7 (0) “ X*{p)) f° r a typical [0/902]* laminate. The conclusions that can be drawn 
from this and other similar plots are as follows: 

1. At low crack densities, microcracking is preferred. At some critical density, the energy release rate for 
delamination will surpass the energy release rate for microcracking and delaminations will be expected 
to initiate at the tips of the microcracks. Once delaminations begin, microcracking will cease and the 
delaminations will grow. 

2. The critical density depends on the laminate structure. For a fixed number of 0° plies, the more 90° 
plies there are, the lower will be the crack density required to get delamination. For a fixed ratio of 0° 
to 90° plies, the more plies there are, the lower will be the crack density required to get delamination. 

3. If the relative toughnesses for microcracking and delamination are different because of material prop- 
erties of because of different deformation modes, the critical density for initiation delamination will 
change. For example, as the delamination fracture toughness gets lower than the microcracking fracture 
toughness, the critical crack density to initiate delamination will get lower. In the limit of relatively low 
delamination fracture toughness, delamination will begin after the first microcrack. 

When delamination is not through-the-thickness, the problem is no longer two-dimensional and requires 
a three-dimensional analysis. Although a complete three-dimensional variational mechanics analysis seems 
intractable, it appears possible to construct a quasi-three-dimensional analysis based on the two-dimensional 
analysis and a lumped spring model of a laminate containing an arbitrary profile of delamination growth. 
This problem will be considered in a future publication. 
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Figure 2: The microcrack density as a function of applied load in Hercules 3501-6/AS4 carbon fiber/epoxy 
cross-ply laminates. The symbols are experimental data points and the smooth lines are best fits using 
G mc of 240 J/m 2 . 
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Figure 4: The x direction tensile stress in the 90° ply in a [90 2 /0], carbon fiber/epoxy laminate (<7rx ). 
This ply has microcracks located at p — ±3. 
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Figure 5: The microcrack density as a function of applied load in a Hercules 3501-6/AS4 carbon fiber/epoxy 
[902/0], laminate. The symbols are experimental data points and the smooth line is calculated using 
G mc of 240 J/m 2 . 
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